import NBody from "NBodyFramework"

s = NBodySystem.Read(argv[0])
cosmo = Cosmology.Model.SCDM
Analysis.VirialRadius.CalculateR200(s, cosmo.Density(cosmo.Redshift(s.Time)))
